The kinetics of escaping of Brownian particles from a potential well for different space 

dimensionality. The effect of external force. 
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The kinetics of two (2D) and three (3D) dimensional diffusion-assisted escaping of Brownian 
particles from a potential well in the presence of an external force is analyzed in detail. The kinetics 
is studied within the two-state model (TSM) proposed for processes in the absence of external force. 
The generalized variant of this model, taking into account the force effect, is proposed which is 
shown to be quite accurate for some shapes of the well both for 2D and 3D processes. Within 
the generalized TSM simple expressions for the well depopulation kinetics and, in particular, for 
the escape rate are obtained. The effect of the force (F) is shown to manifest itself in the escape 
rate dependence on the only parameter <p = Fa/(2kbT), where a is the Onsager radius of the 
attractive part of the well U(r), defined by the relation |f/(o)| ~ kf,T. The limiting behavior of this 
dependence in the cases of weak and strong force is studied in detail both in 2D and 3D processes. 
Some applications of obtained results to the analysis of experiments are briefly discussed. 
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PACS numbers: 82.20.Db, 82.20.Mj, 61.20.Lc 
I. INTRODUCTION 

The effect of external force on mechanisms and kinetic 
properties of condensed phase diffusion-assisted reaction 
processes is considered in a large number of works both 
experimentally and theoretically^ 2 -! 3 -! 4 - The active inter- 
est of scientists to this phenomenon results from its great 
practical importance. 

One of the most important systems, in which the ex- 
ternal force effect is investigated very thoroughly, is re- 
combining geminate ion pairs, undergoing relative diffu- 
sion in the external electric fieldi 2 -! 3 -! 5 -! 6 - Most of theoretical 
studies analyze the kinetics of the recombination process 
within the simplest model, which reduces the problem to 
solving the Smoluchowski equation for probability dis- 
tribution function (PDF) of particles diffusing in a pure 
Coulomb potential (with an external force) and reacting 
with the rate highly localized at short distances^! 7 -* 8 -! 9 - 
Even in this most simple formulation the problem can, 
in general, be solved only numerically, though detailed 
analytical analysis of some simple variant of the prob- 
lem have also been made* 2 -! 7 - for example, within the pre- 
scribed diffusion approximation* 1 ^ 

Recent advances in time resolved investigations of 
charge transfer and escaping processes in fast gemi- 
nate reactions^ 1 - and, in particular, geminate recombi- 
nation of ion pairs in non-polar and moderately polar 
liquid a 11 ' 12 ! 13 ' 14 ! 15 ' 16 ! 17 ' 18 inspire further development of 
theoretical methods of the analysis of the considered 
problem. The main challenge of the theoretical studies 
consists in the correct description of the manifestation of 
specific features of the interparticle interaction (in real 
liquids) in the reaction kinetics in a tractable form sim- 
ple enough to be suitable for applications. 

In the majority of above mentioned theoretical works 
no specific features of the form of the interaction poten- 
tial for the probe (Brownian) particles at short distances 



(of order of molecular size) have been taken into account. 
In the condensed phase, however, the distance depen- 
dence of the potential at short interparticle distances r 
can be strongly modified by interaction of particles under 
study with those of the medium.— This modified interac- 
tion is usually characterized by the so called mean force 
potential (MFP), which in a physically reasonable form 
incorporates the medium effect and, in particular, dis- 
creteness of the medium at short distances. The interac- 
tion with the medium particles is known to result in the 
wavy behavior of the MFP at short distances. Moreover, 
in some systems the medium effect results in the well- 
type shape of the MFP at short distances (see Fig. 1) 
with the a markedly high barrier at distances r of order 
of the distance of closest approach d. This effect is found, 
for example, in the case of ion pairs in polar liquids^ ' 19 ' 20 

Concerning the applicability of well-type approxima- 
tions for the real MFPs, it is also worth mentioning the 
additional reason: from mathematical and kinetic points 
of view any attractive potential can be considered as well- 
shaped in the absence of (or low) reactivity of particles 
at r ~ d. The only difference of this type of wells from 
those shown in Fig. 1 is in their urge-like shape at r ~ d. 

The well-type shape of the MFP (with the reaction bar- 
rier at r ~ d) results in the formation of the quasiequi- 
librium state within the well, which can be considered 
as a ca ge. In the absence of external force the kinet- 
ics of diffusion-assisted depopulation of the initially pop- 
ulated cage state is analyzed in detail in a number of 
papers i 21 ' 22 ' 23 ' 24 In the limit of deep well depth the prob- 
lem is shown to be accurately described with the two- 
state model (TSM), i.e. the model of two kinetically cou- 
pled states: the quasiequilibrium localized state within 
the well and the free diffusion state outside the welli^ 3 . 

The TSM enables one to obtain the well depopulation 
kinetics in a relatively simple analytical form. This kinet- 
ics, determined by the monomolecular reactive passing 
over the barrier at r ~ d and escaping from well (cage), 
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II. FORMULATION OF THE PROBLEM 
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FIG. 1: The picture of the interaction potential u(r) for two 
models of its shape: narrow (dashed) and wide (full) well; 
a n and a w are the Onsager radii for these models [defined by 
u(a„) = 1, {y = n,w)], d is the distances of closes approach, 
and rb is the coordinate of the bottom. 



appears to be non-exponential, in general. In the limit of 
deep well, however, the deviation from the exponential 
kinetics is shown to be fairly smal h 22 ' 23 

In this paper we generalize the TSM to describe the 
effect of the external force on the well depopulation ki- 
netics. The effect is assumed to result only from the force 
induced change of the rate of escaping from the well and 
will concentrate on the discussion of the manifestation of 
the force in the escaping kinetics. 

Within the generalized TSM we derive simple formu- 
las for diffusion assisted well depopulation (or escaping) 
kinetics in two (2D) three (3D) dimensional spaces for 
different models of the well shape. Moreover, taking into 
account that in the limit of deep well (which is of main 
interest of the discussion) the kinetics is close to expo- 
nential, special attention is paid to the analysis of the 
escape rate. The analysis shows that the specific fea- 
tures of the force effect on the escape rate depends on 
the well shape. This effect can be characterized by the 
rate dependence on the only parameter. In the cases of 
weak and strong force the limiting analytical expressions 
for this dependence are obtained and briefly discussed. 

Some possible applications of obtained results are also 
discussed. It is shown that the kinetic formulas are fairly 
useful for the analysis of different experiments. As ex- 
amples we considered recent time resolved experiments 
on liquid phase chemical reaction a 11 ' 12 ' 13 ' 14 , transient 
photocurrents arising from dissociation of exciplexes in 
solution o 16 i 17 ' 18 i 25 , and the kinetics of colloidal particle 
trapping into and escaping from the optical force induced 
potential well (tweezers) 26 ' 27 . 



We consider the dynamics of the Brownian particle un- 
dergoing (2D or 3D) diffusive motion in a spherically 
symmetric MFP well U(r) = U(r) centered at r = 0, 
where r is the vector of position of the particle. The de- 
pendence of U(r) on the distance r = |r| is schematically 
shown in Fig.l. The MFP will be characterized by three 
parameters: the distance a (called hereafter the Onsager 
radius), at which U(a) « fc^i 1 , which is defined as: 

a= Ur° drr ~ 2eU(r) ) ^ 

[in Fig. 1 a n and a w denote the Onsager radii, corre- 
sponding to two models of the well shape, discussed be- 
low (in Sec. V)], the distance of closest approach d, and 
the radius r b of the bottom of the well whose energy is 
U{r b ) = -U b = -{k B T)u b . At r - d the MFP U(r) is 
assumed to be of the shape of a barrier, diffusive passing 
over which models the reaction within the well. 

The main purpose of the work is to analyze the ki- 
netics of diffusion-assisted escaping from the well in the 
presence of the external force — F, i.e. escaping from 
the well of the potential Uf — U(r) + (F • r). For 
dehniteness the force is assumed to be directed along 
the axis z: F = (0, 0,F). The analysis can conve- 
niently be made in spherical coordinates in which r = 
(r sin 9 cos 0, r sin 9 sin 0, r cos 9) . 

The diffusive space-time evolution of the Brownian 
particle is described by the PDF p(r, t) = p(r, 9, <f>\t). In 
general, in spherical coordinates the PDF depends on all 
three variables r, 9, and (f>. However, in the considered 
case of isotropic diffusion and isotropic initial condition 
p(r, t = 0) = Pi(r) the PDF p{r,t) is independent of the 
azimuthal angle <p so that p(r,t) = p(r, 9\t). In our anal- 
ysis we will assume that particles are created within the 
well at initial distance T{ ~ r b : 

p(r,t = 0)= Pl (r)=N t 5(r~r % ), (2.2) 

where N, = (27rr t )" 1 and N t = (47rr 2 )" 1 for 2D and 3D 
processes, respectively. 

The PDF p(r, 9\t) satisfies the Smoluchowski equation 

p = V r {D(r){V rP + pV T Uf)}, (2.3) 

where V r is the gradient operator, 

u t (r) = u(r) + (f • r), (2.4) 

with u(r) = U(r)/(k B T) and f = F/(k B T), is the di- 
mensionless MFP, and D is the diffusion coefficient for 
the particle, which in our analysis is assumed to be in- 
dependent of r [though some possible effects of D(r)- 
dependence can also be studied]. 

Note that the reaction kinetics for pairs of interacting 
Brownian particles, say a and b, is described by the equa- 
tion similar to eq. (|2.3p with r = r a — rb and parameters 
expressed in terms of those for separate particles. 2,6 
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In the absence of force the kinetics of escaping from the 
spherically symmetric short range potential well u(r) is 
analyzed earlier . 21 ' 22 ' 23 ' 24 Here we extend the approach 
applied in these works to describe the effect of external 
force. This approach is based on the approximate solu- 
tion of eq. (|2.3j) in the limit of deep well, in which this 
equation can be solved analytically by expansion in the 
small parameter r c /r e <C 1, where 

r c ~ (a - df/D and r e - T c e u » (2.5) 

are the time of equilibration within the well and the time 
of escaping from the well, respectively 22 ' 23 

Analysis of this solution shows^ 3 - that in the lowest 
order in the parameter r c /r e the Smoluchowski approxi- 
mation (|2.3j) is equivalent to the model of two kinetically 
coupled states: the state within the well and the state of 
free diffusion outside the well. 



III. TWO-STATE MODEL. 3D PROCESSES. 

Originally, the TSM was proposed to treat the kinetics 
of diffusion-assiistcd escaping from the well in the case 
F = 0, in which the PDF p(r, t) depends only the dis- 
tance r = |r|. In the presence offeree [f = F/(fcsT) ^ 0], 
however, the potential Uf (r) = u(r) + (f • r) in eq. (|2.3[) . 
is anisotropic, which results in the dependence of the so- 
lution p(r,t) on particle polar angle 8: p{r,t) = p(r,8\t). 

It is important to note that in the case f ^ the TSM 
is also valid for a variety of shapes of the potential well 
u(r) though some additional analysis of the correspond- 
ing validity criteria are certainly required. 

A. Kinetic equations 

In accordance with rigorous consideratio n 22 ' 23 , in the 
lowest order in r c /r e ~ e~ Ua <C 1 the escaping kinet- 
ics can be described within the model of two kinetically 
coupled states: highly localized state within the well and 
free diffusion state outside the well. The evolution of the 
state within the well (d < r < a) is determined by the 
well population 

n(0\t)=4TT [ drr 2 p(r,8\t), (3.1) 

Jd 

while the evolution of the state outside the well (r > a) 
is governed by the PDF c(r,0|t). Within the TSM the 
effect of the force f shows itself in the dependence of 
kinetic parameters on 8. The form of this dependence is 
determined by the particular variant of the model. 

In general, TSM kinetic equations, describing evolu- 
tion of PDFs n(6\t) and c(r, 9\t) in the presence of an 
external force, can be written as; 22 ' 23 

h = S a K + (8)c(a\t) + [L c - (K-(0) + w r )]n, (3.2a) 
c = L f c+ [S^K-^n- K + (9)c]6(r -a), (3.2b) 



where S a — Ana 2 , 

L f = £>V r (V r + f) (3.3) 

is the operator, describing 3D-diffusion outside the well, 
and L c is the Smoluchowski operator in {#}-space which 
controls orientational relaxation of the PDF in the well. 

The essential difference of anisotropic equations from 
isotropic ones consists in the orientation dependence of 
rates, K + (8) and K-(6). In the considered limit r c /r e <C 
1 we will assume the transition rates K± to satisfy the 
relations^ 3 - 

K ± -> oo and K-{8)/K+{8) = K e (0). (3.4) 

Therefore in this limit ^-dependence of rates K± (8) show 
itself in that of the the equilibrium constant: K e (9). The 
form of the function K e (9) is determined by the shape 
of the well. Some model well shapes and corresponding 
K e {8) dependences, as well as applicability of the corre- 
sponding TSMs, are discussed below. 

Equations (|3.2[) should be solved with boundary con- 
ditions 

(V r + /cos0)c| r=a = and c\ r -,oo — + 0, (3.5) 

first of which describes reflection of particles (diffusing in 
the state outside the well) at r — a. The initial condition 
is assumed to be isotropic and given by eq. (|2.2p . 

In what follows it will be convenient to represent func- 
tions n{8\t) and c(r,8\t) in the form of vectors |n(t)) and 
\c(r,t)}, whose components are obtained by expansion 
of these functions in the orthonornmal basis of properly 
normalized Legendre polynomials (spherical functions) 28 

10 = (Z + ±)P/ (cos 6>), (Z|= / d6 sin 8P t (cos 8)... (3.6) 

Jo 

with I = 0,1,- •• : 

OO OO 

\n)=J2ni\l) and |c)=^ C/ |0, (3.7) 

1=0 1=0 

where for any vector |x(#)), (x = n i c )i its components 
Xi are defined by 

Xi = (l\x) = I d8 S m8P l (co S 8) X (8). (3.8) 
Jo 

In term of this vector representation the initial condi- 
tion can conveniently be written in the form, explicitly 
displaying its independence of orientation: 

|p l ) = (2 7 rr 2 )- 1 |0) ( 5(r-r l ). (3.9) 

As for the initial condition, it is worth noting, in ad- 
dition, that in the most realistic limit of orientational 
relaxation within the well much faster than the escaping 
from the well the escaping kinetics is insensitive to the 
orientational dependence of the initial condition. 
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B. Escaping kinetics 

Equations (|3 . 2[) can be solved by the method applied 
in the absence of force but with the use of expansion of 
n(0\t) and c(r,9\t) in spherical functions |Y/), i.e. vector 
representation |n(t)) and |c(r|t)) [see eq. (|3.7p ]. The 
solution yields for the Laplace transform 

/>oo 

|n(e)) = / dte- et \n(t)) : (3.10) 
Jo 

|S(e)) = [e + w r - L c + VF e (e)]" 1 |nO- (3-11) 
In this expression 

w r = (D/Z w )( f drr- 2 e u ^) . (3.12) 

is the rate of first order reaction in the well in which 

Z w = j drr 2 e' u{r] (3.13) 

J d<r<a 

is the partition function for the well u(r) = [/(r)/(fcsT) 
[in the absence of force (/ = 0)] ) 22 ' 23 

W e (e) = G- 1 (a,a\e)K e (3.14) 

is the operator of the generalized e-dependent escape rate 
(see below), in which the equilibrium constant K e is the 
operator, indicating its dependence on the angle 9, and 
G(a, a\e) is the evolution operator for diffusive motion 
outside the well (evaluated at r — ri — a): 

G(a,a\e) = (a|(e - 

= e-" cose (a|(e- A / )- 1 |a)e vcos9 , (3.15) 

where ip = fa/2 and 

k } = D(L r + r- 2 L e -\f 2 ), (3.16) 

is the auxiliary operator in which 

U = — V e (sin0V e ) =-y2l(l + l)\l){l\. (3.17) 

31110 U 

and 

L r = Dr- 2 V r (r 2 V r ) (3.18) 

are the operators of free orientational and radial diffu- 
sion, respectively. 

Noteworthy is that for the initial condition (|3.9[) |rij) = 
(1/2010). 

In what follows we will restrict ourselves to the analy- 
sis of the well depopulation kinetics n(t), whose Laplace 
transform is given by 

%(e) = 2tt [ d9 sm0n(6,e) = 2vr(0|n(e)) 
Jo 

= (0\[e + w r - Lc + Weie)}- 1 ^). (3.19) 



C. General formulas 

Formulas (|3.15I(|3.17p allow us to evaluate the operator 
G _1 (a, a\e) in analytical form and, therefore, analyze the 
behavior of no(e) relatively easily. 

In the evaluation it is worth taking into account the 
useful relation which simplifies the differential operator 
in the radial space: 

(a\(e-A f r 1 \a) = (a\(e-\ f y 1 \a), (3.20) 

where 

\ f =D(V 2 r +r- 2 L e ~y 2 ). (3.21) 

The evolution operator (a|(e— A/) _1 |a) can be obtained 
in analytical form^S with the use of two linearly indepen- 
dent operator solutions ip-(r) and ip+(r) of equation 

(e - \ f )i)± = (3.22) 

in which the operator Lg is treated as a parameter. These 
solutions satisfy two boundary conditions corresponding 
to those given in eq. (|3.5p [after change of variable c(r) = 

e -(/rcos0)/2^ r )] 

(V P + \f^-\r=a = and ^ + \ r ^oo -> 0, (3.23) 
where 

oo 

u>= \l) (l\ cos 6\l') {I' \ (3.24) 

l,l'=0 

is the matrix representation of the function cos 9. The 
matrix elements (Z|cos0|i') = (Z|Pi(cos ff)\l') are evalu- 
ated analytically 2 ^ though the corresponding formulas 
will not be needed in our further analysis. 

Both solutions ip+( r ) an d V 7 — ( r ) are expressed in terms 
of matrices of Bessel functions K l+ i (x) and I l+ i (x)$& 

oo 

X(r) = VfJ2\ l ) X i+h&er/a){l\, (X = I,K): (3.25) 

?/3+(r) = K{r), -ij)_(r) = t{r) + K(r)k, (3.26) 
where tp = fa/2, 

Ve =^l + e/e f , e f = Df 2 /4=(D/a 2 )<p 2 , (3.27) 

and k is the matrix determined by the boundary condi- 
tion at r = a [see eq. (|3.23p ]: 

k = [V r K{r) - qk(r)]- l [qi{r) - V r /(r)]| r=Q , (3.28) 
in which 

q = a _1 (l - tpoj) = a _1 (l - </?cos0). (3.29) 
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It is worth noting that the matrices K and I do not 
commute with Cj and, therefore, the order of matrices in 
the products of the matrices in expressions (|3.26p - (|3.28p 
is important. As a result of these special commutation 
properties of the matrices, the matrix solutions ip+(r) 
and V>_(r) do not commute either. 

The representation of the evolution operator (r|(e — 
•V ) ~~ 1 1 r i ) m terms of non-commuting solutions ip+ ( r ) and 
(r) is proposed and thoroughly discussed in ref . [29] . 
This representation generalizes the well known one for 
scalar solutions ip+(r) and ijj-(r). In general, the pro- 
posed representation is fairly cumbersome. In the par- 
ticular case of solutions given by eq. (|3.25|) . however, it 
reduces to a more simple one: 

(r\(e-X f )- 1 \r i )=g(r,r i )+K(r)kk(r i )W^ (3.30) 
In this formula 

W m = D[\7 r I(r)K(r) - \7 r k(r)I(r)} = D (3.31) 
is the Wronskian of two solutions and 
g(r,n) = [K(r)I(n)6(r - n) 

+I{r)K(n)B{n - r)]W~} (3.32) 

is the evolution operator for k = in which 0h(x) is the 
Heaviside step function. 

The validity of the expression (|3.30p can be verified by 
direct substitution to equation inhomogeneous variant of 
cq. (|3.22|) with delta-function in the right hand side. 

For the particular case r = r, = a formula (|3.30|) yields 

(a\(e-X f )- 1 \a)=D- 1 [q + q K (e)]- 1 , (3.33) 

where 

oo 

q K {e) = -V r k(r)/K(r)\ r = a = ]T \l)q Ki m\ (3-34) 

1=0 

with 

q Kl (e) = a' 1 [I + ip e Ki_i (<Pe)/K l+ i (<p c )] . (3.35) 

Substituting the expression (|3.33p into eq. (|3.20p and 
then into eqs. (|3.16p and (|3.14p we obtain formula 

G- l {a,a\e) = D[q + e' vcos0 q K (e)e lpcos6 }. (3.36) 

For our further analysis of the escaping kinetics no (t) 
we need to specify of the operator L c describing orien- 
tational relaxation in the well. Naturally it should be of 
the Smoluchowski-like form: 

L c = ^(sin^Velsm^Ve + V e u)], (3.37) 

where D c ~ D/r% is the orientational diffusion coeffi- 
cient u(8) is the effective orientational potential which is 
determined by the shape of the well (see below). 

Moreover, in the considered limit of large well depth it 
is quite natural to assume that orientational relaxation 
is much faster than well depopulation. 



D. Fast orientational relaxation in the well 

The fast orientational relaxation limit implies that 
D c 3> t~ . This relation means that after some time 
~ t c = D^ 1 of orientational relaxation (of the initial pop- 
ulation in the well) the vector of well population |n(i)) 
remains close to the equilibrium one |^ e ) during the pro- 
cess: 

|n(t)> w n(t)|* e ), (3-38) 

where 

|* e ) = Z^e-^W, Z g = [ d0sm8e- a » (e \ (3.39) 

Jo 

Note that within bra-ket notation the adjoint vector (ip e \ 
coincides with (Yq \ and is given by formula 

(* e | = (0| = / dO sin^..., (3.40) 
Jo 

which can be confirmed by the relation {^ e \L c — di- 
rectly following from the definition of L c [see eq. (|3.37p ]. 
With the use of this formula one can easily find that \ip e ) 
satisfies the normalization condition (^ e \^ e ) = 1. 

In what follows we will restrict ourselves to the analysis 
of the escaping kinetics just in this limit of fast orienta- 
tional relaxation. 

For fast orientational relaxation the splitting SL C of 
eigenvalues of the operator L c (SL C ~ D c ) is much larger 
than || Well ~tu°. In such a case in the lowest order in the 
parameter w®/D c C 1 we can significantly simplify the 
general expression for n (e) [eq. (|3.19p ] and thus for the 
inverse average lifetime wo — Tq 1 = Uq (0) as follows: 

n(e) = [e + m r +w e (e)] _1 and wo = w r + w e , (3.41) 

where 

w e (e) = (* e |We(e)|* e ) and w e = w e (0). (3.42) 

Equation (|3.4ip presents the main result of the work 
for the kinetics of the well depopulation in the limit of 
fast orientational relaxation. 

The depopulation kinetics predicted by formulas (|3.4ip 
and (|3.42p essentially depends on the shape of the po- 
tential well which determines the orientational poten- 
tial u(0) in the Smoluchowski-type operator L c [see eq. 
(|3.37p ] and, therefore, the equilibrium state |W e ). In our 
work we will consider two realistic models of the well 
shape in which simple analytical expressions well depop- 
ulation kinetics can be obtained. 



E. Application of results 

In this section we will analyze the specific features of 
mean escape rate w e for two variants of the well shape: 



6 



1) Narrow well shape (shown in Fig. 1 by dashed line 
with a — a n ), for which a — d <C d and the time of 
equilibration within the well r c ~ (a — d) 2 / D <C a 2 / D <C 
T e , where r e is the time of escaping from the well [see eq. 

2) The wide well shape (full line in Fig. 1 with a = 
a w ), corresponding to a small distance of closest approach 
ti«« (in which of the main interest is the region fd <C 1 
while fa> 1). This shape is schematically shown in Fig. 
1 by full line with a = a w . 

The analysis will be made in the above-discussed limit 
of fast orientational relaxation in the well with the use 
of eq. (|3.41|) for the well depopulation kinetics. In our 
study, in addition to this formula we will also use the rep- 
resentation for e ±¥,cos9 in terms of expansion in spherical 
functions P;(cos( 



.30 



3 ±<p cos t 



- £(±l)^ + ±)/ ;+ ^)P ; (cos#). (3.43) 

in 



1=0 



1. Narrow-well shape 

In the case of narrow well, when a — d <C d, the well is 
of the shape of attractive well layer near the distance of 
closest approach d. In this limit within the wide region 
force strengths / < 1/ (a — d) we can neglect the effect 
of the force on the radial shape of the well and take into 
consideration only the dependence of well depth u/O) on 
the orientation angle 6: 



Ui(9) s» Uf (rft, 9) ~ ujj + fa cos 9 



(3.44) 



with / = |f| > 0, and the force effect on free diffusion 
in the state outside the well. In eq. (|3.43|) we took into 
account the smallness of the width of the well, a — d <C d, 
which leads to the high accuracy of the relation frb ~ fa- 
It is important to note that the small value of the 
well width and, therefore, fast equilibration of the well 
population in radial direction, ensures the validity of 
the description of the kinetics in terms of the angular 
coordinate dependent well population n(9\t) introduced 
above. Noteworthy is also that the negligible force af- 
fected change of the well shape results in the absence of 
the dependence of the detailed balance relation and the 
equilibrium constant K e (9) on the angle 9. In such a case 
K e {9) is given by the relatio n 22 ^ 3 



K e {9) = K° e = a 2 /Z u 



(3.45) 



in which the effective partition function Z w is given by eq. 
(|3.13[) . i.e. is controlled by the shape of the potential u(r) 
without external force, despite possible strong force effect 
on the energy of the bottom predicted by eq. (|3.44p . This 
is because the external force leads to the identical change 
of both the bottom energy Ub{9) ~M|, + fa cos 9 and the 
energy of the free diffusion state at r = a: U{(a,9) w 
fa cos 9. 



The potential Ub{9) determines the kinetics of orienta- 
tional relaxation of the population in the well, which is 
described by the Smoluchowski operator (|3.37[) with 

u(9) = Ub{9) — Ub — 2ipcos9, where ip = fa/2. (3.46) 

In this case the equilibrium state within the well is writ- 
ten as 

I*., . Wltt ZoW . !i£MM. (3 . 47) 



Substitution of formulas (|3.4T|) and (|3.43[) into the ex- 
pression (|3.42p yields for the function w e (e) 

w e (e)/w° en =Q(e) = \ + ^coth(2^) + S{ip, e). (3.48) 

In this formula 22 - 23 



w° e = w e {<p = 0, e = 0) = Da/Z u 



(3.49) 



is the escape rate in the absence of a force and 
2tt 



ipZ ((p) 



X^ + ^+ifeOafog. (3.50) 



1=0 



where 

Qiife) =aq Kl (e) = I + y e K x _ i ((p e )/K l+ i (<p e ). (3.51) 



with (p e = ipy/1 + e/ef [eq. ([3~27| ], 

The function w e (e) — w®Q(e) is fairly complicated so 
that, in general, the depopulation n(t) can hardly be ob- 
tained in analytical form. In the considered limit of deep 
well, however, the main e dependent contribution ~ (p e 
comes from the term with I = of the sum in eq. (|3.50[) . 
The e-dependence of other terms with I > 1, which are 
of higher order in ip e (~ y> 2 ), can be neglected by taking 
ip t w tp e= o = tp. In so doing one gets 



S(tp,e) fa S(ip) + cp 1 ta,nh(tp)((p e - tp), 



(3.52) 



where: 



%>)=%>,0) = 



and 



2-Kif' 



.\ oo 



1=0 



w e (e) ^w e + w v [(l + e/e f ) 1/2 - 1] (3.54) 
with w e = w e (e = 0) = w°Q(0) [see eq. (j3~4^1) ]. 

Cf = (D/a 2 )tp 2 , and w v — io e tanhy. (3.55) 
The inverse Laplace transformation of n(e) yields 2 ^ 



m = [ 

2m J_ 



V (3.56) 



D+0 l + e + 7(e/ + e) 1/2 
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where £/ = ef/w, 



w 



w v , and 



7 



w® tanh ip 



7 e with 7 e = 



(3.57) 



The kinetics (|3.56p can be expressed in terms of error 
functions^ In the absence of force formula (|3.56[) . natu- 
rally, reduces to that obtained earlier for ip = ^ 



a =o(t) = 



1 

2tt7 



ioo+0 exp[ £ ( Wo t)] 

l + e + 706 1 / 2 ' 



(3.58) 



ioo+0 



where Wq — w((p = 0) = w r -\- and 70 = 7^=0 — 

«/ U ;o) 1/2 7e. 

Here we are not going to discuss general properties of 
kinetics n{t) but restrict ourselves to the qualitative anal- 
ysis of its asymptotic behavior at short and long times in 
the limit 7 e <C 1 , corresponding to the case of deep poten- 
tial well, and in the most realistic case of relatively weak 
force, in which Ef < 1, i.e. D/w < 4// 2 . In this limit at 
relatively short times t < w^ 1 ln(l/j v ) ~ w^ 1 \n(l/j ) 



the kinetics is exponential; 22 ^ 3 n(f) 



-wot 



(w 



w r + w e ). In the opposite limit t w~ ln(l/7 v ) the de- 
pendence n(t) is non-exponential i 22 ' 23 nit) ~ i~ 3 / 2 e _e/t . 
As expected, in the absence of force (e/ =0) the kinetics 
at long times becomes of inverse type power, which is de- 
termined by long trajectories of particles diffusing freely 
outside the welL 22 i 23 Note that for weak force ip <C 1 the 
parameter e/ is small: e/ ~ ip 2 , and can be neglected. 
This means that in the weak force limit with high accu- 
racy the kinetics fiit) coincides with that for ip = in 
which w = wq. Detailed analysis of specific features of 
the kinetics predicted by formulas (|3.56|) and (|3.57|) is 
made in ref. [24]. 

It is easily seen that that in the considered limit of 
deep well the the non-exponential inverse power tail of 
the kinetics no(t) is small. Moreover the effect of this 
tail tends to zero as the value of the force is increased. In 
such a case the effect of force ip on the escaping kinetics is 
reduced to that on the mean escape rate w e (<p) = w e (e = 
0,V>) 



Ml* 



n(v) = h+<pcoth(2<p) + S(ip). (3.59) 



Formula (|3.59|) shows that the force effect on the rate 
is characterized by the only parameter p = fa/2. The 
numerical calculated universal function Q n (tp) which de- 
scribes this effect is displayed in Fig. 2a. In addition, 
some limiting specific features of the behavior of Q n (p) 
can be revealed with simple analytical expressions. 

a. Weak force limit. In the limit of weak external 
force, when (p = fa/2 <C 1, in two lowest orders in p 
the mean escape rate w e can be estimated using only the 
first term (with I = 0) in the sum S(tp) in the expression 

for Q n (p) =We(p)/w° en : 



Q n {p) = Q^{p)~ 1 + ^ = 1 + fa/2. 



(3.60) 



b. Strong force limit. In the opposite case <p = fa/2 
1, denoted as the strong external force limit, the analysis 
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FIG. 2: The force (tp — fa/2) dependence of dimensionless 
escape rates Q v (tp) = w e (cp)/w ev (p) (y = n,w) [for narrow 
(n) and wide (w) wells], calculated with exact eqs. (|3.59[) 
and (|3.66[) (full lines) and interpolation formulas (|3.62p and 
(|3.71[) (triangles). Shown also is the dependence of the di- 
mensionless difference of rates 8Q w (ip) [see eqs. (|5.ip and 
(|5.2[) ] calculated with the exact formula (full line) and the 
interpolation expression (|3.71|l (squares). 



of the dependence w e {ip) with the use of eq. (I3.59|) is 
more complicated. However, the corresponding limiting 
dependence can easily be derived taking into account that 
in the limit (/>> 1 the escaping process is, actually, one- 
dimensional, for which the Q(y>)-factor is given by 2 ^ 

Q n (p) = Q s n (p) ~ 2<p = fa. (3.61) 

In deriving eq. (|3.61[) we used the expression for one- 
dimensional escape rati 
2.(1) - f ,/„-«('■) . 

.7 Jd<r<a U ' e 

partition function. 

Note that this relation is also confirmed by the semi- 
quantitative estimation with the use of eq. (|3.53p by 
truncating the sum at I = I* ~ ^JTp 3> 1 and then sum- 
ming up the terms which mainly contribute to S(ip). 

c. Interpolation formula. The limiting results obtained 
above for ip -C 1 and ip> 1 can be combined into a simple 
algebraic interpolation formula 

Qn\f) 



g± w x = Df/Z£\ in which 
« Z w /a 2 is the one-dimensional 



Qn(f) « Wn VP) = J--t-<^V^ _ e ') 

which reproduces function Q n ((p), numerically evaluated 
using eqs. (|3.53|) and (|3.59[) . with accuracy ~ 3% (see 
Fig. 2a). 



l + p{2-e' lp ) 



(3.62) 
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2. Wide-well shape 

Another form of the well shape, in which analysis of 
the escape rate w e can be made analytically, corresponds 
to the small distance of closest approach, or large a, for 
which d ~ Tf, -C a. In this case in a fairly wide region of 
relatively strong force / < 1 /d, 1 /r& the escape kinetics 
is fairly accurately described by the TSM (|3.2b[) . 

It is important to note that the inequality fr b < 1 en- 
sures quite high accuracy of the approximation neglect- 
ing the effect of force on the well shape in the region near 
the bottom. In this approximation, the quasiequilibrium 
population distribution within the well is isotropic: 

|*e) - |0). (3.63) 

This, in turn, means that the partition function Z w is 
independent of the angle 9 and is given by eq. (|3 . 1 3|) . 

The effect of force, however, manifests itself in the 
anisotropy of the activation energy of escaping u a {9): 

u a (0) « ut(9,a) w u b + 2<pcos0, (<p = fa/2), (3.64) 

which, in turn, leads to the anisotropy of the detailed 
balance relation, i.e. the anisotropy of the equilibrium 
constant 

K e {6) = K°e~ 2lpcos6 , (3.65) 

where K® is the isotropic equilibrium constant in the 
absence of external force given by eq. (j3.64|) . 

Formula (13. 65|) calls for some additional comments es- 
pecially concerning its applicability. The fact is that 
the value of K e {ff) at each particular 9 is determined 
assuming local quasiequilibrium of the population out- 
side and inside the well in the region close to r = a at 
this 9. In general, it is difficult to justify the existence of 
the quasiequilibrium in the considered limit, unlike the 
limit of narrow well discussed above. This is because for 
ip = fa/2 < 1 the time of passing over the escaping bar- 
rier width §b ~ min{a, (the width of the region of 
transition from the inner part of the well to the outer 
one), Tb ~ d 2 /D is comparable with the time of reorien- 
tation r c ~ a 2 /D. It is worth noting, however, that the 
accuracy of quasiequilibrium assumption becomes better 
with increasing / since the for ip — fa/2 ^> 1 the width 
6b <C a and, correspondingly, 1% -C t c . 

The above-mentioned arguments lead us to the con- 
clusion that in the considered limit of small radius of the 
well bottom the TSM with ^-dependent equilibrium con- 
stant K e {9) gives quite reasonable interpolation formula 
for the kinetics of the escaping process and, in particular, 
for the escape rate w e , which correctly describes both the 
limit of weak and strong external force. Further analysis 
(see below) will confirm this statement. 

Formula for the escaping kinetics can straightforwardly 
be derived with the use of general formulas (|3.4ip , (|3.42p , 
and some results obtained above in the limit of narrow 
potential well. The fact is that, in the mathematical 



form, the average of any operator multiplied by angu- 
lar dependent equilibrium constant (|3.65[) [of type of 
eq. (|3.42p ] over the isotropic equilibrium state (|3.63p 
is similar to the average over the equilibrium distribu- 
tion (|3.47p . except for the partition function Zo(ip) [eq. 
(I3.47P ]. which should be replaced by Z (p — » 0) = 2 cor- 
responding to the isotropic distribution. These simple al- 
gebraic manipulations result in the following expression 
for w e (tp) 

u> e (¥>)/<. (<P) = QtofrO = \Z^p>)e~ 2 *Q n (p>) (3.66) 

where Zo(ip) and Q n (ip) are determined in eqs. (|3.47p 
and (I3.59p . respectively, and 

0)= We V (3.67) 

is the escape rate in the absence of the external force 
but with the activation energy it*, corresponding to the 
orientation 9 = tt (most favorable for escaping): 

u* a = u a {9 = n)=u b - 2<p. (3.68) 

As in the case of narrow potential well the dependence 
of w e on the force / is expressed in terms of that on the 
only parameter ip. The characteristic function Q w (f), 
which determines the pre-exponential factor in the acti- 
vation type dependence of w e (<p), is displayed in Fig. 2b. 
The numerical results show that Q w (>p) monotonically 
decreases (with increasing ip) from Q w = 1 at <p = Q to 
Qw — 1/2 at if — > oo. This behavior is markedly differ- 
ent from that of Q n (<p) although the specific features of 
Qu> (^-dependence are essentially determined by those of 
Q n {<p)- Some of features of the function Q w (cp), for ex- 
ample saturation at <p — > oo, looking unexpected at first 
sight, can be understood by simple analysis (see below). 

a. Weak force limit. In the weak force limit the behav- 
ior of Q w (ip) at (f = fa/2 <§; 1 differs form that obtained 
above for Qnif) (i-e. for narrow potential well): Q w (ip) 
decreases with increasing ip, so that at ip -C 1 

QS(V)=Q»(^«1)«1-V- (3-69) 

Such a behavior of Q w (p) results from using the in- 
dependent normalizing rate w®^ ~ e 2<p (instead of W® = 
u>e) in the definition of Q w (ip). 

b. Strong force limit. In the opposite limit ip = 
fa/2 ^S> 1 the force strongly affects the average escape 
rate w e , first of all, because it significantly changes the ac- 
tivation energy of the rate w e . As for Q w (ip), which char- 
acterizes the pre-exponential factor of the corresponding 
Arrenius-type expression for w e , at <p ^> 1 it monotoni- 
cally decreases approaching the asymptotic value 1/2. 

The obtained Q w (^-independence at ip — > oo can eas- 
ily be understood by taking into account that, according 
to formula (|3.65p . in the case of wide well for strong exter- 
nal forces the flux of escaping particles is highly localized 
in a small region of orientations 59 = tt — 9 < 1/yftp 
1. The escape rate is determined by the total flux J e 
through this region of size s e ~ {&9) 2 ~ ip -1 . In the 
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strong force limit tp ^> 1 the flux J e ~ tp, as it follows 
from eq. (|3.6ip , so that Q w {tp) ~ s e J e ~ const. The ex- 
act estimation can be obtained just by substitution of the 
corresponding limiting expression (|3.61[) into eq. (|3.66[) : 

Q s w {tp) = \Z^)e-^Q s n ^)\ v>>1 = | (3.70) 

Similarly to the narrow well limit, in the case of wide 
well for large tp the escape rate is determined by the 
quasi-one-dimensional flux of escaping particles. The 
mechanism of formation of the one-dimensional flux is, 
however, somewhat different in both cases: for narrow 
wells the transition to the one dimensional regime re- 
sults from high localization of the well population in the 
small region at 9 ~ n, while for wide wells this transition 
is caused by strong localization of favorable transition 
rates in this region. 

c. Interpolation formula. A simple interpolation ex- 
pression for Q w {tp) can be derived, for example, with the 
use of similar formula for Q n {tp) presented in eq. (|3.62p : 

Qw(<p) « = \Z Q (<p)e-*r& n (<p). (3.71) 

Quite satisfactory accuracy (~ 3%) of this formula is 
demonstrated in Fig. 2b. 

IV. ESCAPING KINETICS IN 2D CASE 

In this section we briefly discuss the kinetics of force 
affected two dimensional escaping of Brownian particles 
from the well. 

The formal analysis of the 2D problem is very similar to 
that in the 3D cases made above. The differences reduce 
to a few changes of analytical forms of dependences on 
the parameters of the model. 

First, the 2D TSM equations are similar to 3D eqs. 
(j3.2p but with 2D variants of parameters and opera- 
tors: S a — 2na, Lf — Z?V r (V r + f) is the opera- 
tor describing 2D diffusive motion outside the well, and 
L c = D c \7g(\7 s + \7 e u a {9)) is the operator of orienta- 
tional (circular) diffusion in the circular cage in 2D space. 
The analytical form of the equilibrium constant K e {9) is 
also somewhat different from those discussed above in 3D 
case (see below) i 22 i 23 

Similar to 3D case the orientation dependent popula- 
tion of the well n(9\t) and PDF of particles outside the 
well c(r, 9\t) are conveniently represented as vectors by 
expanding them in basis of vectors 

\l) =Nicos(W), (l\ = [ d9cos(W)..., (/ > 0), (4.1) 

J — 7T 

where N t >i = l/n and N = 1/(2tt). 

The solution of (2D) TSM equations for vectors |n(t)) 
and |c(r|t)) by the method applied in the 3D case yields 
for n {e) — J* 7r d9n{9,e) = (0|n(e)) the expression 
(I3.19[) . In the fast orientational relaxation limit this ex- 
pression, naturally, reduces to formula (|3.41[) . In both 



these equations the operator W e (e) is given by eq. (|3.14|) 
in which 

G~ 1 (a, a\e) = D\e- V cos e q K (e)e* 3 cos 9 - § / cos 9] , (4.2) 

where q K {e) = YmIo with 

q Kl (e)=a- 1 [l + p e K l _ 1 (ip e )/K l {ip e )]. (4.3) 

With the use of thus obtained general expressions and 
the relation^ 

e±V cose = h{if) +2 ^;~ i J 1 ( ¥ ,)cos(W) (4.4) 

one can derive relatively simple formulas for the depopu- 
lation kinetics no(t) within two models of the well shape 
discussed in Sec. HID: the models of narrow and wide 
well. 

It is worth noting that the kinetics of 2D diffusion- 
assisted reactions are essentially different from that of 3D 
processes because 2D diffusion is recurrent while 3D one 
is transient. In the particular case of diffusion-assisted 
escaping from the well the recurrent nature of 2D diffu- 
sion manifests itself in strongly non-exponential escaping 
kinetics even in the case of deep well. Moreover in the 
absence of external force w e = w e (e = 0) = (4' e |W e (e = 
0)| x I , e ) = 0, i.e. the probability of escaping from the well 
is zero and escaping is always reversible (which is quite 
natural for recurrent process). 

This, in turn, means that unlike 3D escaping processes, 
in the absence of force, the efficiency of 2D ones cannot 
be characterized by the only parameter. The external 
force, however, results in the irreversible escaping flux 
which shows itself in w e {tp ^ 0) ^0 and nearly exponen- 
tial escaping kinetics approaching the exponential with 
the increase of force. Below we will restrict ourselves to 
the evaluation the mean escape rate w e {p) describing the 
force effect on the escape process, which, as we have al- 
ready mentioned above, is getting closer to exponential 
with increasing the force (ip). 

A. Narrow- well shape 

In the 2D case of narrow wells a — d <C d, which is 
defined similarly to the 3D case (see Sec. III.E.l.), in 
the large region of forces / the force effect reduces to the 
change of the quasiequilibrium distribution within the 
well. The corresponding 2D quasiequilibrium state \^ e ) 
can be written as 

|* e ) ^Io 1 {2p)e~ 2ificos6 . (4.5) 

With the use of this expression and formulas, derived 
above, we get 

w e {<p)/w* n = Q n {p) = ph{2p)/I a {2p) + S p {cp), (4.6) 
where 

w* en = D/Z&\ with = [ drre-< r \ (4.7) 

J d<r<a 
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is the auxiliary parameter of dimensionality of rate and 

%)=^ - 1 (2^)^J l / i 2 Mif ( _ 1 ( V )^ 1 M, (4.8) 

with Co = 1 and 0>i = 2. 

The numerically calculated dependence w e ((p) (|4.6[) is 
displayed in Fig. 3. As expected w e (<p — ► 0) — > 0. In the 
opposite limit (/;> 1 we get the dependence 

w e (<p > l)/< n » 2 V = fa (4.9) 

similar to that obtained in the 3D case and corresponding 
to the ID escaping mechanism (see Fig. 3). To clarify 

this statement note that for narrow wells zffl rs aZ^\ 
where Zw — fd <r<a dr e~ u ^ is the ID partition func- 

(2) 

tion, and therefore w e (ip^§> 1) ps 2<pw* e = fa{D/Z w ') m 
Df/ZL 1] . 



B. Wide- well shape 

The opposite limit of (2D) wide wells (a — d ^> d), 
defined by analogy to the 3D case (Sec. III.E.2.), can 
be treated with the use of formulas obtained in the 
narrow-well case. To do that one should take into ac- 
count, similarly to the 3D case, that in the wide- well 
limit |* e ) = |0) and K ( 9) ~ e - 2¥,cose . This means (as it 
is mentioned in Sec. III.E.2) that both limits the aver- 
age w e (e) = (\l/ e } |W 7 e (e)|\E r e ) is similar from mathematical 
point of view and differ only of the normalization factor 
1^(2?) 

w e (ip)/w* ew = Q w {<p) = Q n (<p)I Q (2cp)e- 2lfi , (4.10) 
where 

<=<e 2 ^. (4.11) 

The dependence w e ((p) (|4. 10[) is shown in Fig. 3. In 
agreement with the above statement in the limit of weak 
force w e ((p ~ > 0) — > 0. In the opposite limit of strong 
force ip 3> 1, however, one gets w e (p) ~ yftp as it follows 
from the analysis of the numerically calculated function 
w e ((f) presented in Fig. 3. The dependence w e (p ^> 
1) ~ J~<p> can easily be understood by estimations similar 
to those, clarifying the u) e ((/^-independence of ip at ip 3> 
1 (see Sec. III.E.2). The fact is that for p > 1 the 
density escaping flux J e ~ ip is localized is localized in the 
small sector of size 86 ~ 1/ ^pp around the most favorable 
direction 9 = tt. Therefore the total escaping flux Q w ~ 
J e S6 ~ JTp. The exact asymptotic dependence w e {p S> 
1) can be obtained with eq. (|4. 10[) and formula (|4.D|) : 



w e (<p > l)/w* w » yfpfa. (4.12) 

High accuracy of this expression at ip > 1 is demonstrated 
in Fig. 3. 



V. DISCUSSION AND APPLICATIONS 

A. General remarks 

This work concerns detailed theoretical study of the ef- 
fect of the external force / = Fj (fcsT) on the kinetics of 
diffusion-assisted depopulation of a deep isotropic poten- 
tial well in 2D and 3D processes. Fairly simple matrix ex- 
pressions for the depopulation kinetics are obtained and 
thoroughly analyzed. 

In our work we have concentrated on the analysis in 
the most physically reasonable limit of fast orientational 
relaxation of the population in the well. In this limit 
the analytical expression for the depopulation kinetics is 
derived which predicts the kinetics to be close to the ex- 
ponential in the wide region of parameters of the model 
both for 2D and 3D processes. The total depopulation 
rate in this case is shown to be a sum of reaction and 
escape rates. In our work we have mainly studied the 
specific features of the escape rate w e whose value ap- 
pears to significantly depend on shape of the well. Sim- 
ple analytical expressions for w e (f) are obtained for two 
limiting types of wells: narrow wells of type of well layer 
at a distance of closest approach d (for which a — d <C d) 
and wide wells with large effective Onsager radius a» d. 

In the case of narrow well the effect of the force on the 
escape rate is fairly strong but shows itself only in the 
preexponential factor of the Arrenius-type dependence of 
the rate, i.e. no strong effect on the activation energy is 
predicted. On the contrary, in the case of wide well (or 
small distance d) the force affects not only preexponential 
factor but the activation energy as well. 

It is worth noting that the effect of an external force on 
the diffusion-assisted processes in the presence of interac- 
tion between particles are studied in a number of works 
(see, for example, refs. [1] and [7]). Especially compre- 
hensively the force effect (electric field effect) is analyzed 
in the case of ion pair recombination reaction, i.e. in the 
case of the Coulomb interaction between particles. 

Unfortunately it is practically impossible to compare 
the results of our analysis with majority of those ob- 
tained earlier, since these works mainly concerned with 
processes in potentials without well at short distances, 
the reactivity is usually assumed to be high. In partic- 
ular, in the case of ion pair recombination processes the 
recombination kinetics is considered to be determined by 
diffusive motion in the pure Coulomb potential^^ It is, 
nevertheless, interesting to note that in the small field 
limit fa <C 1 the force effect on the probability P e (f) 
of escape from the Coulomb potential, found in ref. [5], 
is independent of the initial distance between ions and 
is represented in the form P e (f) ~ Pe(f = 0)(1 + fa/2), 
which is in apparent agreement with the field dependence 
of the escape rate obtained in our work [see eqs. (|3.60p 
and 10311]. 

Noteworthy is also that in the strong force limit the 
escaping process becomes nearly one-dimensional in both 
cases of well shape considered in both for 2D and 3D 
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processes. In this limit the escape rate is determined by 
the ID flux in the small region of favorite orientations 
corresponding to 8 ~ ir. This fact allows one to easily 
improve the considered TSM, in which the effect of the 
force on the location of top of the barrier (assumed to 
be at r = a) is neglected. Moreover, in the strong force 
limit one can also take into account the smoothness of 
the shape of the realistic barrier near the top, which in 
the TSM is actually assumed to be of cusp shape. 

Concluding this short discussion of results we would 
like to note that in this work we restricted ourselves to 
the analysis of the most realistic limit of fast orienta- 
tional relaxation within the well. In reality, however, 
with the use of general formula (I3.19|) one can also de- 
scribe the manifestation of finitcncss of the orientation 
relaxation time. The case, in which the effect of finite- 
ness is largest, of course, corresponds to L c — 0, i.e. the 
absence of orientational relaxation. In this case the angu- 
lar dependence of the equilibrium rate [K e (8)] 7 evidently, 
results in the highly non-exponential well depopulation 
kinetics no(t), which can be approximated by the sum of 
exponentially decreasing (monomolccular) contributions 
with ^-dependent rates, coming from different orienta- 
tions. With the use of obtained formulas there will be no 
difficulties to analyze this case as well, when needed. 

In our further discussion possible applications of for- 
mulas obtained we will consider some types of processes 
recently actively analyzed experimentally. 



plied to the analysis of kinetics of some liquid-phase 
reactions ,2iL§i2 Concerning the results obtained above 
for the case / ^ 0, note that the most convenient for 
experimental analysis is not the force dependent inverse 
mean lifetime u>o(v?) [see eq. f|3 .42[) ] . but the difference 
wo((p) — I0o(0) = w e ((p) — ffljj, which is independent of 
the rate w r of reaction within the well (assumed to be 
independent of tp) . The corresponding dimensionless pa- 
rameters 

6Q v (<p) = [m{<p)-M0)]/K,(<P)> (» = n > w )> (5- 1 ) 

are directly related to Q„(y)'- 

SQnfr) = Qn{<p) - I, SQwiv) = Qw{<p) ~ e~ 2(p - (5-2) 

The behavior of SQ n (p) is, clearly, similar to that of 
Qniv) except for evident displacement along ordinate 
axis. As for i5Q u ,((/?)-dependence, shown in Fig. 3, its 
form is essentially different from that of Q w (ip): at ip — > 
the function 8Q w (tp) ~ <p is similar to SQ n (ip), while 
5Q w (tp — > oo) = 1/2. Moreover 5Q w (ip) has a maximum 
(though not very pronounced) at ip — ip m w 2.0. 

It is of great interest to compare these theoretical pre- 
dictions with experimental results of type of those given 
in refs. [9-13] but in the presence of electric field. 

C. Transient photocurrents. 



B. Condensed phase reaction kinetics. 

The obtained formulas are quite suitable for the anal- 
ysis of diffusion-assisted condensed phase geminate reac- 
tions. The effect of external force on escaping and reac- 
tion yields is of special interest in the particular case of 
reaction ion-pair recombination reactions in polar solids 
and liquids, in which the external force can be realized 
by applying electric field. There are a number of exper- 
imental time resolved spectroscopic studies of reactions 
with the participation of ions . n i 12 ' 13 i 14 ' 15 

Traditionally the results of such investigations are an- 
alyzed with the use of model calculations in which in- 
teraction is assumed to be pure Coulomb [.]. It is worth 
noting, however, that in polar media the medium affected 
interaction, which can be described by the MFP, strongly 
deviates from the from the Coulomb one [.]: unlike the 
Coulomb potential the MFP oscillates at short distances 
of order of molecular size [.]. At distances close to that of 
the first coordinate shell the MFP has a most deep well, 
whose depth can be much larger than fc^T for dielectric 
constants e > 20. In this case the proposed approxima- 
tion of the realistic interaction with the use of well-type 
potential is much more accurate than pure Coulomb po- 
tential approach. 

Some of earlier results of the proposed TSM, concern- 
ing the kinetics of processes in the absence of exter- 
nal force (/ = 0), have already been successfully ap- 



Recently considerable interest has been attracted to 
experimental investigations of the recombination kinet- 
ics of ion pairs in non-polar and moderately polar 
liquids with time resolved measurements of transient 
photocurren t 16 i 17 ' 18 i 25 . The photocurrent is found to be 
fairly sensitive to the spatial evolution of photoinduced 
ion pairs and, in particular, to the recombination kinet- 
ics. 

The proposed theory can be very fruitful for the in- 
terpretation and description of experiments on transient 
photocurrents. The fact is that the TSM considered 
above enables one to quite accurately describe the spatial 
evolution of ion pairs without solving the Smoluchowski 
equations which are rather complicated even for restric- 
tive models of the interparticle interaction. In this short 
discussion we will outline and illustrate the possibilities 
of the proposed method. 

Recall that the TSM is valid in the limit of the time t c 
of the PDF relaxation in the well much shorter than the 
average lifetime fo = w^ 1 in the well [according to eq. 
(|3.4ip Wq 1 ~ (tUg) -1 , lUjT 1 ]. Noteworthy is also that the 
TSM describes the kinetics at relatively long times t > 
t c . In reality, however, the method can straightforwardly 
be extended to also treat the PDF relaxation at short 
times t < t c (i.e. relaxation within the well) just because 
of assumed essential difference between thermalization 
time t c and To- This difference leads, in fact, to the 
independence of the processes of intrawell relaxation and 
well depopulation, which can be described separately. 
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The initial stage of relaxation of the ion-pair PDF 
within the well results in the transient photocurrent 
J e (t) at times t < r c . The relaxation kinetics depends 
on the well shape. Below, for simplicity we will con- 
sider the case of narrow well and assume that the initial 
PDF in the well is isotropic. For narrow wells the in- 
trawell PDF relaxation can, in turn, be separated into 
two stages: radial and orientational relaxation. The first 
stage is the radial relaxation whose characteristic time 
is t w = A^/D, where A w — a — d < d, is shorter than 
the time r c = D~ x ps r 2 jD of the second stage of ori- 
entational relaxation to the anisotropic quasiequilibrium 
ion-pair PDF in the well of the potential itf (r) (|2.4[) . 

Together with the final stage, described by the TSM, 
the kinetics of the transient current relaxation can be 
represented as follows: 

1. At t r~j t w (the stage of radial relaxation in the well) 
the transient current can approximately be described by 
the simple expression 

Je(t)« Jo + (Jo- Jo)e-* /r ™, (5.3) 

where Jo = eDf is the current produced by freely diffus- 
ing ion pair [initial current for the isotropic initial PDF 
(12. 2[) ] and Jo§Jo is the final current after radial relax- 
ation [see eq. (I5.4| 

2. At t ~ r c (the stage of orientational relaxation in 
the well) the current J e (i) can be evaluated using the 
relation J e = p(t), where p(t) — e(rf(t)) is the the aver- 
age electric dipole moment of ion pairs within the well, 
in which where r/ is the projection of radius vector of 
the ion pair on the direction of the field: rt = (r • f/)/. 
For the case of narrow well the derivative of the aver- 
age dipole moment p(t) is easily evaluated in the limit 
of weak external force / in the lowest order in / (in the 
linear response approach): 

J e (t) » Joe' 20 "* with J = | J , (5.4) 

where Jo is defined in eq. ()5.3|) . and D c ss D/r 2 (r& is 
the radius of the well bottom). 

3. At t > t c (the stage escaping from the well) the 
TSM predicts the following expression for the current: 

J e (t) pa Jon e (t) with n e (t) = n s {t) — n(t). (5.5) 

Here 

n s (t) = 1 — w r / dr n{r) (5-6) 
Jo 

is the ion-pair survival probability. 

Taking into account different time scales of these three 
kinetic stages one can describe the current relaxation ki- 
netics by the combined expression 

je(t) = Je(t)/J « |(2 + e -*/r w)e -2J3 c t + ^ (g ?) 

This simple expression enables one to quantitatively 
describe the specific features of the current relaxation 



kinetics in a wide region of times. In deriving eq. (|5.7p we 
have assumed the well to be narrow. As applied to the ion 
pair recombination and transient current relaxation this 
assumption is quite appropriate in the limit of relatively 
polar solvent with dielectric constant e s > 15, in which 
the ion-pair MFP u(r) is known to be of the shape of 
deep narrow well at short distances r ~ rii 19 i 20 

Formula (|5.7|) . however, turns out to be of quite rea- 
sonable accuracy even in the case of diffusion-assisted 
processes in the pure Coulomb potential u{r) — a/r with 
nearly reflective (low reactivity) boundary condition at 
r = d, whose well (at r > d), at first sight, can hardly be 
treated as narrow for a»d. The fairly good accuracy of 
eq. (I5.7[) in this case results from the cusp-like shape of 
the well at r ~ d, where the major of the well population 
is localized (see below). 

To demonstrate the accuracy of eq. (|5.7[) we will com- 
pare its predictions with recent calculations of the current 
relaxation kinetics based on numerical solution of the 
Smoluchowski equation^ These calculations have been 
made for u{r) = a/r, with a = 63.3 A, and for small 
force (p — 0.05, which only very weakly affects the escap- 
ing rate w e . The other parameters of the model are; 2 ^ 
the coefficient D = 3.1- 10~ 5 cm 2 /s of relative diffusion of 
ions, the contact distance d = 9 A, and the reactivity ay, 
which determines the reactive flux at a contact distance 
r = d: D[W r p + (V r u + f cos9)p]\ r=d = p(<r r /4ird 2 )\ r=d , 
and whose value a r = 0.01 • (AwDa) corresponds to the 
weak reactivity limit. 

Recall that the discussed TSM is applicable in the limit 
of small 7e = (c^w^/D) 1 / 2 < 1 [see eq. ([3~F7|) ]. However, 
for the parameters of the model chosen above 7 e ps 0.9, 
i.e. the TSM is not expected to reproduce the exact 
numerical results very accurately. Nevertheless, even in 
this case the accuracy of this model appears to be quite 
good (see below). 

Noteworthy is also that in the considered model assum- 
ing not very deep well of cusp-like shape the quasiequilib- 
rium PDF in the well is somewhat displaced to r > d. In 
this case, with reasonable accuracy one can put rj = f = 
fjjjdrre-^W/Zy, ~ 19 A so that D c « 8 • lO^" 1 [eq. 
(|5.4p ]. For this value of r& one can also approximately 
estimate t~ 1 : t" 1 ps D/(r b - d) 2 ps 3.1 • H^s -1 . 

For the chosen parameters of the model the survival 
probability nf = n s {t — ► oo) w 0.085. 25 Taking into 
account that the TSM exactly predicts the asymptotic 
(at t — > oo) behavior of the escaping kinetics we can 
write the relation nf = w e /uio = w e /(w r +w e ) — 0.085. 

Figure 3 displays the comparison of the time dependent 
survival probability n s (t) and the transient current J e (t) 
with those calculated by numerical solution of the Smolu- 
chowski equation^ For convenience of the comparison 
with numerical results in Fig. 3 we use the dimensionlcss 
time t = Dt/a 2 . The comparison shows good accuracy 
of the TSM prediction for the kinetics n s (t), which is de- 
termined by the evolution of the system at long times 
r > 0.1. As for the time dependence of the transient 
current J e (r), it is also quite accurately reproduced at 
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FIG. 3: The comparison of time (r = tD/a 2 ) dependences 
of the ion-pair survival probability n s (r) and the normalized 
transient current j e (r) — J e (r)/Jo, calculated analytically 
with eqs. (|5.6[l and (|5.7[l (full lines) and numerically 2 ^ (dashed 
lines). The parameters used are presented in Sec. V.C. 

long times, when the contribution of n s (i)-kinetics [rep- 
resented by the last term in eq. (|5.7[) ] is dominated. Some 
(not very strong) difference between analytical and nu- 
merical behavior at shorter times r < 0.1 is a result of 
above-mentioned approximations applied in deriving eq. 
(|5.7p for J e (t) at first two short time stages of the current 
relaxation. 



D. Optical tweezers. 

The important problem which have recently attracted 
much attention, and in study of which the TSM can be 
very suitable is the kinetics of trapping of colloidal parti- 
cles by optical tweezers, i.e. by the optical force induced 
potential wells (arising due to gradients of laser beam 
intensity) j 26 i 27 Detailed experimental investigations show 
that the tweezers potential well u(r) is highly localized 



(the size is about ) and is expected to be highly 
localized.— Of course, in general the well is not spher- 
ically symmetric. Moreover, recently the axial (non- 
potential) component of the trapping force is found^ 2 - 
However, to a good accuracy, for description of trapping 
and escaping kinetics one can neglect these non-trivial 
contributions to the trapping force and use the isotropic 
approximation for the tweezers well. 

The proposed theory based on the TSM is very useful 
for the analysis of the kinetics of trapping into and es- 
caping from tweezers potential well. Application of the 
TSM significantly simplifies the problem of description of 
the kinetics thus allowing for the analysis of more com- 
plicated effects of a large number potential wells. 33 

Of special interest is the effect of well motion on the 
trapping/escaping kinetics discussed in a number of pa- 
pers concerning possible applications of tweezers.— This 
effect is known to reduce to that of the external force dis- 
cussed above. The fact is that the well motion induces 
the force F acting on a particle. In the frame of refer- 
ence, moving with the well, the force is proportional to 
the well velocity v: F„ = /i~ x v, where /x = Z?/(fcsT) is 
the mobility of the particle, i.e f„ = I? _1 v. Note that 
this relation holds in the case of time dependent velocity 
v(t) as well. 

Thus the problem of the analysis of the effect of well 
motion is equivalent to that thoroughly discussed above 
with the external force f = f v (and without reactivity, i.e 
with w r — 0) . In our consideration we have assumed that 
f is independent of time. The case of time dependent f 
[for example, because of time dependent velocity v(t)] 
is, in general, much more complicated and can hardly be 
analyzed analytically. However, in a quite realistic case 
of relatively slowly changing force, for which the charac- 
teristic changing time r/ > 1/wo, one can treat the force 
effect adiabatically evaluating the nearly exponential es- 
caping kinetics with the use of formulas derived above for 
static /, in which the escaping rate w e ((p) is replaced by 
the corresponding time dependent expression w e ((p(t)). 
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